Interacting Fibonacci anyons in a Rydberg gas 
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A defining property of particles is their behavior under exchange. In two dimensions anyons can 
exist which, opposed to fermions and bosons, gain arbitrary relative phase factors p] or even undergo 
a change of their type. In the latter case one speaks of non-Abelian anyons - a particularly simple 
and aesthetic example of which are Fibonacci anyons [2]. They have been studied in the context of 
fractional quantum Hall physics where they occur as quasiparticles in the k = 3 Read-Rezayi state 
[3], which is conjectured to describe a fractional quantum Hall state at filling fraction v — 12/5 
[4] . Here we show that the physics of interacting Fibonacci anyons [5] can be studied with strongly 
interacting Rydberg atoms in a lattice, when due to the dipole blockade [S] the simultaneous laser 
excitation of adjacent atoms is forbidden. The Hilbert space maps then directly on the fusion 
space of Fibonacci anyons and a proper tuning of the laser parameters renders the system into an 
interacting topological liquid of non-Abelian anyons. We discuss the low-energy properties of this 
system and show how to experimentally measure anyonic observables. 



There is great interest in studying many-body systems 
of interacting anyons as they often exhibit exotic quan- 
tum phases [THS]. A further motivation is that the 
exchange of anyons - the braiding - permits the imple- 
mentation of robust protocols for quantum information 
processing [TDHT^] . Physically, anyons emerge as quasi 
particles on the ground state of an interacting many- 
body system and recently there has been put much ef- 
fort in implementing and exploring anyonic models with 
cold atoms [T5HTo] , polar molecules [H] and trapped ions 
|17j . A model which has been much studied - in partic- 
ular in the context of quantum information processing 
- is Kitaev's toric code [13]. Excitations on the ground 
state are Abelian anyons, which merely acquire a phase 
when braided. A more complex and exotic scenario is en- 
countered in the non-Abelian case, i.e. when anyons can 
undergo an actual change of their type under braiding. 

A particularly simple example of non-Abelian anyons 
are Fibonacci anyons which occur in two types: A trivial 
particle referred to as 1 and a non-trivial particle denoted 
by r. A convenient way to define an anyonic system is 
through fusion rules for the particle types (see e.g. Refs. 
[2 [12] ) , which for Fibonacci anyons read 

1x1 = 1, 1xt = t, t X 1 = r, rxr = l + r. 

The first three rules are a consequence of the trivial na- 
ture of a 1-anyon. The fourth rule states that two anyons 
of type r can fuse such that the result is either the triv- 
ial particle or a r-anyon. This can be thought of as be- 
ing analogous to the merging of two spin 1/2 particles 
{\ ® \ = © 1), which can yield either a singlet (total 
spin 0) or a triplet (total spin 1). 

For a set of N anyons of type r it can be shown that 
a basis of the underlying Hilbert space is spanned by all 
possible fusion paths (2J, i.e. the number of ways in which 
these anyons can be successively fused. For large N, this 



number grows as ip with ip = (1 + \/5)/2 being the 
golden ratio. An example of a fusion path is illustrated 
in Fig. [T^i and it is constructed as follows: Starting from 
the left the first r-anyon fuses with the r that is set 
by the left boundary condition. There are two possible 
fusion outcomes, 1 or t. In this particular realization 
we choose the outcome 1, which is written into the first 
box. This trivial particle then fuses with the second r- 
anyon. The fusion rules dictate that the outcome must 
be a r anyon (second box). Continuing this procedure by 
obeying the fusion rules results in a possible fusion path. 
It is important to note that since a trivial particle fusing 
with a r always results in another r the occurrence of 
two consecutive l's is excluded. 

In the following we will show in detail that a Ryd- 
berg lattice gas constitutes an analog quantum simula- 
tor for Fibonacci anyons. The link between these two 
systems is the aforementioned exclusion principle. While 
for anyons this is a consequence of the underlying mathe- 
matical rules the exclusion in a Rydberg system is phys- 
ically rooted in the dipole blockade [5j which prevents 
the simultaneous excitation of neighboring atoms to Ry- 
dberg states. The blockade originates from large elec- 
trostatic energy shifts between atoms in Rydberg states 
and has recently been demonstrated experimentally for 
atoms trapped in individual optical traps [THl HO] ■ Also 
first experiments with Rydberg atoms in a lattice |21j 
have been conducted. We will discuss which anyonic in- 
teractions can be realized in such a system and how any- 
onic degrees of freedom can be measured experimentally. 
We expect that this comparatively simple quantum sim- 
ulator platform for non-Abelian anyons will highlight a 
new route towards the study of exotic forms of quantum 
matter. 

Formally, the Rydberg lattice gas is described by a 
set of atoms with the electronic ground state \g) = \l) 
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FIG. 1: (a) Fusion space. The Hilbert space of a system of 
N Fibonacci anyons of type r is spanned by all possible fusion 
paths. An example of such path is given by the symbols con- 
tained in the boxes. A fusion path translates into configura- 
tions of atoms which are either in the electronic ground state 
or excited to a Rydberg state. The Rydberg state (ground 
state) is identified with the fusion outcome 1 (r). (b) Level 
scheme and interaction energy. The excitation of a Ry- 
dberg atom on site k shifts the energy of Rydberg states of 
atoms in the neighborhood. The strong interaction between 
neighboring atoms (here k and k + 1) excludes their simulta- 
neous excitation to Rydberg states and effectively constrains 
the Hilbert space of the atomic system to the fusion space of 
Fibonacci r-anyons. 



that is coupled by a laser with Rabi frequency 51 to a 
high lying Rydberg state denoted as |r) = Excited 
atoms with position labels k and m interact with a van- 
der-Waals potential V] fc _ m | = (Ce/a 6 )/\k — m| 6 where 
C*6 is the interaction dispersion coefficient and a the lat- 
tice spacing. Within the rotating-wave approximation 
the Hamiltonian of an ensemble of N — 1 atoms reads 
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Here A is the detuning of the laser frequency with re- 
spect to the frequency of the atomic transition \g) <+ |r), 
rifc — (1 + o~ k : )/2 is the local number operator and tx^' 2 
are Pauli spin matrices referring to the internal degrees 
of freedom of each atom. Hamiltonian (fTl) has been stud- 
ied in a number of theoretical works [22H25| and recent 
experiments have shown that it accurately reflects the 
physics of laser driven Rydberg gases [2TJ . 

We are here interested in a parameter regime in which 



the nearest neighbor interaction is the largest energy 
scale, i.e. |f2|, |A| <c |Vi|. Fig. [T]d depicts a sketch 
of this situation where we take V\ (C 6 ) [23] negative. 
However, in this work we will discuss both positive as 
well as negative V\. It was shown in Refs. [331 US] that 
when |f2|, |A| <C \V\\ the dynamics of the Rydberg gas is 
confined to a Hilbert space "Hbiockado that is spanned by 
atomic configurations in which adjacent atoms are not si- 
multaneously excited. This Hilbert space is equivalent to 
the fusion space of an ensemble of N Fibonacci anyons of 
type t. To see this, one should think of the atoms as be- 
ing located in between the r-anyons. A fusion path (see 
Fig. [iji for an example) is then encoded in the internal 
state of the atoms: the fc-th atom is in the state |r)(|<?}) if 
the outcome of the fusion of the A:-th anyon with the pre- 
vious fusion outcome is 1(t). The state of each atom can 
thus be interpreted as the combined topological charge 
of the anyons located to its left [5]. 

Having established a formal equivalence between the 
Hilbert spaces of a Rydberg lattice gas and the fusion 
space, we will now turn to the question concerning which 
anyonic interactions are actually realized by the atomic 
system. To this end, it is more convenient to work with 
an explicit projection of Hamiltonian (JlJ on the con- 
strained Hilbert space H blockade- This effective Hamil- 
tonian reads H c s = Hq + -ff v dw + H2 with 
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where = 1 — n&. Here Hq + i? v dw is the actual pro- 
jection of Hamiltonian ([I]) onto "Hbiockadc- The term H2 
contains second order corrections that arise from a non- 
perfect blockade, i.e. due to a finite nearest-neighbor 
interaction. Since the relative strength of H 2 scales with 
(fi/Vi) 2 it can be regarded as a perturbation to Hq. 
Moreover, due to the quickly decaying tail of the van- 
der-Waals interaction also the term -ff v dw can be con- 
sidered perturbatively. We will later discuss the influ- 
ence of these terms on the spectrum, but let us first an- 
alyze the principal Hamiltonian Hq. We start by rewrit- 
ing it as Hq = ^2 k=2 2 hk + H\y where the h k are local 
three body Hamiltonians and H^ defines the boundary 
term. Using the basis (centered around the k-th site) 
{Itlt) , lilt) , I til) , lltl) , llll)} and taking into account 
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that these site-triples overlap, we find that 
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where /i = A + V 2 -/3 + 27, / 2 = A-a-/3 + 37 and a, /3 
and 7 are constants. These constants are arbitrary, show- 
ing that there is a large class of local three-body Hamil- 
tonians whose sum adds up to Hamiltonian H . The 
boundary term reads ifg = ( A — a + 7)711 + (fz — 7)^2 + 
(a-7)n A r_2 + (A-/ 2 -7)n A r_i + rj [crfP 2 + fV-20%_i] 
for open and = hx+hj^-i for closed boundaries (with 
relabeling the indices N — > 1, — > N — 1). 

An anyonic interaction which has been extensively dis- 
cussed in the literature 2 , 5 is the anyonic equivalent of 
the Heisenberg interaction between spins. As mentioned 
earlier two spin 1/2 particles fuse according to the rule 
I ® I = © 1, where and 1 are the total spins (sin- 
glet and triplet) and the Heisenberg interaction, oc Si-S 2 , 
introduces an energy difference between these fusion out- 
comes. According to r x t = 1 + t an anyonic analogue 
thus would assign different energies to the two fusion out- 
comes 1 and t. The construction of this Hamiltonian is 
discussed in Refs. [H[S], but we will briefly summarize 
here how an explicit representation of the interaction in 
fusion space and hence also in "Hbiockade is obtained. 

When considering two r anyons there are different or- 
ders in which anyons can be fused, two of which are 
depicted in Fig. [2^. The upper panel represents a 
"conventional" fusion path according to the rules dis- 
cussed earlier: Here the left anyon of type x merges 
with the first r-anyon. The outcome of this fusion is 
an anyon of type y which fuses with the second r to 
yield an anyon of type z. This is not a convenient ba- 
sis for the discussion of anyonic interactions, as we can- 
not directly read off the fusion outcome of the two r's. 
We thus perform a basis change from the "conventional" 
fusion space | xyz) e {|1t1) , \tt1) , \1tt) , \t1t) , \ttt)} 
to a basis where the fusion outcome (denoted as 
u) of the two r-anyons becomes explicit: \xuz) € 
{|111> , \tt1) , \1tt) , \t\t) , \ttt}}. The change to this 
new basis, as depicted in Fig. [2^,, is conducted by the 
so-called F-matrix [21 [5] 
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With this we are now equipped to construct the Heisen- 
berg interaction. An operator that discriminates the two 
fusion results is a projector on the trivial particle. In 
the basis where the fusion outcome of the two r's is ex- 
plicit (bottom panel of Fig. [2^,), this is a diagonal ma- 
trix with entries one, where u — 1 and zero otherwise, 



Uk = diag(l, 0, 0, 1, 0). Upon transforming Ilk back into 
the "conventional" fusion space we can thus write the 
anyonic Heisenberg Hamiltonian as i? a n = — ^X)*Li ^fc 
with 
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Comparing the coefficients of with those of the three- 
body Hamiltonians hk in Eq. ^ we find that this inter- 
action is in fact naturally present in a system of interact- 
ing Rydberg atoms when 



- J ip 
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-Jitp- 2 -^ 1 ), V 2 = -J<p. (5) 



With this choice of parameters the Rydberg lattice gas 
behaves just like a system of Fibonacci anyons that inter- 
act with a Heisenberg interaction. For J > (V 2 < 0), it 
is energetically favorable for anyons to fuse to the trivial 
particle. In the opposite case, J < (V 2 > 0), the fusion 
to a t- anyon is favored. 

Before turning to the analysis of the spectral properties 
and the influence of imperfections such as a non-perfect 
Rydberg blockade or the van-der-Waals tail let us show 
how one can experimentally measure properties of the 
anyonic system. We have mentioned earlier that the state 
of an atom encodes the combined topological charge of 
the anyons located to its left. This property can hence 
be readily inferred from a site-resolved projective mea- 
surement. Another important experimental quantity - in 
particular in the context of the discussed Heisenberg in- 
teraction - is the result of a fusion of two r-anyons. To 
determine the fusion outcome of the anyons with label k 
and fc + 1 one needs to determine the expectation value of 
the projection operator . This amounts to a measure- 
ment of a three-point correlation function on the atomic 
degrees of freedom. If the outcome is 1(0) we know that 
the particles have fused to an anyon of type l(r). (n^) 
can be inferred from three separate projective measure- 
ments each of which concern only a single site. Such 
site resolved addressing in conjunction with projective 
measurements has been experimentally demonstrated in 
the context of interacting Rydberg atoms [TH1 [20] ■ A 
schematics of the procedure and involved levels is given 
in Fig. [2]}, and a comprehensive list of all measurement 
outcomes together with the corresponding fusion result 
is provided in Fig. [2b. 

Let us finally analyze the ground state and the low- 
energy properties of the model. The ground state is in 
general a strongly correlated many-body quantum state 
whose preparation typically requires an adiabatic pas- 
sage protocol: Here, one starts in a simple initial (prod- 
uct) state with no Rydberg atoms present and slowly 
tunes the laser parameters until one reaches the point 
([5]) in parameter space. This has to be achieved in a 
time much smaller than the radiative lifetime of Rydberg 
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FIG. 2: (a) Changing the order of the fusion path. The F-matrix is performing a basis change in fusion space such 
that the fusion outcome u of the two neighboring r-anyons becomes explicit, (b) Determination of the fusion outcome 
via a measurement of the three-point correlation function (life}. We start by applying a strong laser pulse with Rabi 
frequency Q TB S> | Vz | on a transition from the Rydberg state |r) to a stable hyperfine state \s). This state transfer is necessary 
for switching off the Rydberg-Rydberg interaction. We then measure the state of the atoms on sites k ± 1 by monitoring 
fluorescence from a closed transition, i.e. we shine in a laser that resonantly couples the state \g) to a short lived state \a) that 
decays under the emission of a photon (with rate V) back to \g). If the atom is in state \g) this results in the cyclic emission of 
photons while in the opposite case (atom in |s)) no photons are emitted. The presence/absence of these scattered photons thus 
allows directly to infer the atomic state, (c) List of measurement results and the corresponding fusion outcome. If 
one finds that the atoms located at sites k — 1 and k + 1 are in different internal states we directly conclude from the matrix 
representation Q that (life) = and hence the fusion outcome was r. In turn, if the atoms are in the state | s) fc _ 1 |s) ft+1 , we 
can conclude that (life) = 1 (fusion of the two r's to a trivial anyon). The most involved scenario is encountered when the 
two outer atoms are in the state \g) k _ 1 |<7} fe+1 - In this case we have to measure the expectation value of the operator that 
corresponds to the lower right 2 x 2-block of Fife given in Eq. Q. This block can be written as a rotation of the number operator 
n = (1 + o z )/2: R^nRn. Here R n = e l ' rn ' <T / 2 with a = (a x ,a v ,a z ) conducts a rotation around the axis n = (^~ 1//2 , 0, f^ 1 ) 
by an angle n. Hence, if we rotate the state of the fc-th atom on the Bloch sphere by applying R n and subsequently perform a 
projective measurement we find the fusion outcomes to be l(r) if the atom is in the state |s) ([<?)). 



atoms, which is typically (for principal quantum numbers 
in the range n = 40 ... 60) on the order of 100 /is. The ex- 
perimental feasibility of such protocol has recently been 
shown |27j . Beyond that it is moreover possible to infer 
properties of the ground state such as critical behavior 
from the time evolution and the scaling properties of cer- 
tain quantum mechanical observables as demonstrated in 
Ref. [2£]. 

The ground state phase diagram of the principal 
Hamiltonian H a has been studied in [28] . It has a number 
of interesting features, e.g. for certain parameter choices 
the ground state becomes a simple matrix product state 
[23] and also the excitation gap can be calculated ana- 
lytically [35] • There are furthermore integrable lines in 
parameter space that also include the point ([5]) which 
represents the anyonic Heisenberg interaction [oj. Here 
one can map the model onto the restricted-solid-on-solid 
model {29] and the critical properties of the ground state 
of Hq are inferred from the exact solution. We consider 
the two cases (i) antiferromagnetic interaction between 
anyons (J > 0, V 2 < 0) and (ii) ferromagnetic inter- 
action (J < 0, V% > ), separately. For case (i), the 
low-energy excitation is gapless and the continuum limit 
of the model is described by the minimal conformal field 
theory (CFT) with central charge c = 7/10 30J. For 
case (ii), the low-energy excitation is again gapless but 
the continuum limit is described by another CFT with 



c = 4/5 [30]. The left column of Fig. [3] shows the low- 
energy spectra obtained by exact diagonalization up to 
N = 23 with periodic boundary conditions. For (i), two 
gapless modes at momenta K = and K = tt are clearly 
visible. On the other hand, for (ii), the gapless modes 
are located at K = 0, 27r/3, and 47r/3. These gapless 
modes suggest the existence of the quasi-long-range or- 
der. The difference of modulation period between (i) and 
(ii) can be understood as follows: In the original physical 
system, J > ( J < 0) and hence V 2 < {V 2 > 0) means 
the attractive (repulsive) interaction between atoms in 
Rydberg states. But since there is a hard-core constraint 
caused by the Rydberg blockade, the density wave order 
of period 2 is favored when J > 0, while that of period 3 
is favored when J < 0. 

We finally discuss the effect of the hitherto neglected 
perturbations, i.e., -ff v dw an d H 2 . Fig. [3] shows in the 
right column how the low-energy spectra are modified 
by these perturbations. The spectrum for case (i) with 
-ffvdw an d H 2 suggests that the gap is small and the 
model still remains in the vicinity of the critical point. 
As demonstrated in Ref. [5J, the gaplessness of the an- 
tiferromagnetic anyonic chain is protected by translation 
symmetry and topological (Y) symmetry. Since the first 
and the second terms in H 2 can be absorbed into Hq by 
fine-tuning A and V 2 , we expect that the criticality of 
the model is robust against the longer-range interactions 
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stemming from the tail of the van-der-Waals potential 
and the correction due to the non-perfect blockade. On 
the other hand the spectrum for case (ii) clearly indi- 
cates the existence of a gap. The numerical results do 
not show any evidence of the quasi-degeneracy of the 
ground state caused by the density wave order. We thus 
conclude that the perturbations -ff v dw an d H 2 lead to a 
disordered gapped state that does not break any symme- 
tries. 

In conclusion, we have shown that a Rydberg lattice 
gas rather naturally constitutes an analog quantum sim- 
ulator platform for Fibonacci anyons with tuneable inter- 
actions. Anyonic observables can be monitored by mea- 
suring atomic correlation functions which opens up the 
possibility for probing correlated quantum states of non- 
Abelian anyons in current experiments. 



FIG. 3: Low energy excitation spectra. Energy spectra 
for Ho (left column) and H c g (right column) with positive and 
negative J (referred to as case (i) and case (ii) in the text). 
The spectra have been shifted such that the lowest energy 
eigenvalues are zero. The dashed lines sketch low-energy dis- 
persion relations and serve as a guide to the eye. The different 
symbols (green circles, red squares, blue diamonds) indicate 
different system sizes N: top panels (19,21,23), bottom panels 
(16,19,22). 
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